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Excitation of Surface Waves on Conducting, Stratified, 
Dielectric-Clad, and Corrugated Surfaces 

James R. Wait 

An expression for the field of an electric dipole located over a flat surface with a speci- 
fied surface impedance Z is derived from the formal integral solution by a modified saddle- 
point method. Using the value of Z appropriate for a homogeneous conducting ground, 
the general expressions reduce to those given by Norton. In this case the phase of Z lies 
between and 45°. When the phase exceeds 45°, as it may for a stratified ground, the radi- 
ated wave of the dipole becomes partially trapped to the interface. This effect is most 
pronounced for an inductive surface where the phase of Z is 90°, in which case the energy 
of the wave is confined within a small distance from the surface. Such inductive surfaces 
are a metallic plane with a thin dielectric film or a corrugated surface. 

This unifying treatment provides a link between the surface waves of Zenneck, Sommer- 
feld, Norton, and Goubau, and indicates that the phase angle of Z controls the extent to 
which these waves may exist for a dipole excitation. 

1. Introduction 

The propagation of radio waves along the surface of the ground has been discussed from a 
theoretical standpoint for many years. As long ago as 1907, Zenneck [1] l showed that a wave, 
which was a solution of Maxwell's equations, traveled without change of pattern over a flat 
surface bounding two homogeneous media of different conductivity and dielectric constants. 
When the ripper medium is air and the lower medium is a homogeneous dissipative ground, the 
wave is characterized by a phase velocity greater than that of light and a small attenuation in 
the direction along the interface. Furthermore, this Zenneck surface wave, as it has been 
called, is highly attenuated with height above the surface. In 1909 Sommerfeld [2] solved the 
problem of a vertical dipole over a homogeneous ground (half-space). In an attempt to explain 
the physical nature of his solution, he divided the expression for the field into a "space wave" 
and a "surface wave." Both parts, according to Sommerfeld, arc necessary in order to satisfy 
Maxwell's equations and the appropriate boundary conditions. The surface-wave part varied 
inversely as the square root of the range, and it was identified as the radial counterpart of the 
plane Zenneck surface wave. 

For many years it was believed that this Sommerfeld surface wave was the predominant 
component of the field radiated from a vertical antenna over a finitely conducting ground. It 
was not until 1935 that an error in sign in Sommerfeld's 1909 paper 2 was pointed out by Norton 
[3], which also partly accounted for the unusual calculated field-strength curves of Rolf [4]. 
At about this time there was a series of papers by Weyl, Van der Pol, Niessen, Wise, and Norton, 
deriving more accurate representations for the field of the dipole. A discussion of this later 
work has been given by Bouwkamp [5]. Norton [6, 7], in particular, has developed his formula 
for the field components to a stage where numerical results can readily be obtained. It is now 
generally accepted that the Sommerfeld surface wave (or the radial Zenneck wave) does not 
bear any similarity to the total field of a vertical dipole over a homogeneous conducting (Mirth. 
In fact, the field excited by a dipole varies as 1/d, where d is the distance for low frequencies, and 
it varies as 1/d 2 for high frequencies. Norton has suggested that the field in air of a dipole 
over a homogeneous ground be expressed as a sum of three components: A direct ray (or pri- 
mary influence), a reflected ray which is to be modified by an appropriate Fresnel reflection 
coefficient, and a correction term. Norton has described the first and second components as 
the space wave; the third or correction term, the surface wave. This seems to be a logical step, 

1 Figures in brackets indicate the literature references at the end of this paper. 

2 It has been suggested by B. M. Fannin (Final Report, Pt. Ill, Investigations of air to air and air to ground experimental data, Elec. Eng. 
Dept., Cornell University, 1951) that Summerfeld did not actually make an "error in sign" in his 1909 paper. Fannin's contention is that Rolf and 
other workers simply misinterpreted Sommerfeld's formula for the attenuation factor. This can hardly be so, as Sommerfeld himself in 1910 (The 
propagation in wireless telegraphy, Jahrb. d. Drahtl. Tel. 4, 157, 1910) used the incorrect formula and its resulting series expansion. 
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although, taken separately, the space and surface waves of Norton are not solutions of Max- 
well's equations. On the other hand, his "space wave" is the contribution that would be derived 
on the basis of geometrical optics, and his "surface wave" is the correction from wave theory. 
This latter term will be called the "Norton surface wave" as distinct from the Zenneck and 
Sommerfeld surface waves, and the trapped surface waves discussed below. 

It was pointed out by Sommerfeld [8] nearly 60 years ago that a straight cylindrical con- 
ductor of finite conductivity can act as a guide for electromagnetic waves. Some 50 years later 
Goubau [9, 10] demonstrated that such a cylindrical surface can be launched with reasonable 
efficiency from a coaxial line whose outer surface is terminated in a conical horn. The improve- 
ment in the transmission properties by coating the wire with a thin dielectric film has been dis- 
cussed in detail by Goubau [10]. The plane counterpart of the Goubau-Sommerfeld cylindrical 
surface wave is obtained when a flat metallic surface is coated with a dielectric film. Atwood [1 1] 
has discussed the nature of the surface waves that can exist in a structure of this type. When the 
film thickness is small compared to the wavelength, the single propagating mode has a phase 
velocity slower than that of light and is attenuated rapidly above the surface in the air. Such 
surfaces have been called inductive because the normal surface impedance, looking into the 
surface, is almost purely imaginary for a low-loss dielectric on a well-conducting base. A similar 
type of surface wave can exist over a corrugated surface [12, 13, 14] which is inductive if the 
periodicity and depth of the corrugation are small compared to the wavelength. 

The excitation of surface waves on dielectric-clad plane conducting surfaces has been 
discussed by Tai [15] for a line current source, and Brick [16] for a dipole source. Correspond- 
ing treatments for corrugated surfaces have been given by Cullen [17] for a line magnetic or 
infinite slot source, and Barlow and Fernando [18] for a vertical electric dipole source. In the 
case of dipole excitation, the field varies predominantly as the inverse square root of the dis- 
tance along both the dielectric-clad and the corrugated surfaces. Surface waves of this type, 
which have a phase velocity less than that of light, can be called " trapped surface waves" 
because the}" carry most of their energy within a small distance from the interface [19, 20]. 

It is the purpose of the present paper to study the general problem of a vertical dipole 
over a flat surface with a specified surface impedance, Z. A solution is obtained for the field 
in terms of the error function of complex argument. It is shown that, using an appropriate 
value for Z, the expressions for the Norton space and surface wave for a homogeneous conduct- 
ing ground are obtained as a special case [7]. By allowing Z to be purely imaginary and using 
asymptotic expansions for the error function, the field is shown to be of a trapped-surface- 
wave type [14, 19]. Another special case is for a two-layer conducting ground [21, 22] which 
can exhibit some of the features of trapped surface waves. Thus, a general connection between 
these various forms of seemingly unrelated surface waves is developed. 

2. Basic Formulation 

A vertical electric dipole of length ds and current / is located at a height h over a fiat 
ground plane which exhibits the property of surface impedance. Here the ratio of the tangen- 
tial electric and magnetic fields is specified on the surface of the ground plane and is denoted 
by Z. Choosing cylindrical coordinates, the dipole is located at (0,0, h) and the ground plane 
at z=0. As will be seen, the Hertz vector at (p,<M) need have only a z component, n, to 
satisfy the boundary condition ; consequently, the fields are given by 3 

d 2 IT / b 2 \ 

*>=§£' *=°' M^) 11 (1) 

/3 2 dn 
fffl—0, H <it =- — -=-—j H z =0, 

9 9 inoo dp 

where /x=47rX10~ 7 and ^=2x/free space wavelength. The resultant Hertz vector, II, can now 

3 A time factor exp (iut) is suppressed throughout. 
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be represented as a sum of a primary and secondary pari in the manner 

n=n p +n r , (2) 



where 

Ids 



n^rzrce-^i/r, withc=^- and ri =[(z-h) 2 +p 2 ]^ (3) 



and II r is to be regarded as the influence of the ground plane. In the most elementary case 
of a perfectly conducting ground plane (i. e., Z=0), II r would be given by the simple relation 

U r =ce-^^/r 2 with r 2 =[0+A) 2 +p 2 ]^. (4) 

To find n r (or II) for the more general case, it is necessary to express the primary field in 
terms of an integral representation [2] (sometimes called Sommerf eld's integral), as follows: 

r x Jo " u 

where u= (X 2 — /3 2 ) 1 * and J (\p) is a Bessel function of order zero and argument Xp. Due to 
the presence of the branch point of u, the contour of the integral is indented upward slightly 
at \=p. Physically, this integral can be regarded as a means of representing the primary field 
as a spectrum of plane waves whose angles of incidence, 0, on the ground plane are related to 
the integration variable by A=/3 sin 6. Consequently, the angles include both real and com- 
plex values, ranging from to 7r/2 along the real axis and from 7r/2 to 7r/2+i<» along a line 
parallel to the imaginary axis. 

The form of eq (5) suggests writing n r in the following form: 



■HT* 



CK P )e- uiz+h) R(\)\d\ 



u 

where E(\) is now a reflection coefficient for plane wave incident on the ground plane at angle 
of incidence 0=sin -1 X//3. Consequently, 

/i(X)= ^+^A' (7) 

where 

Z/3 Z 



A = 



/xco 120tt 



In the limiting case of a perfect conductor (Z=0) it can be seen that jf2(\) = l, and the integral 
is of the form of a Sommerfeld integral and II r is again given by eq (4). 
The Hertzian potential II can then be written in the form 



n=c — -+- 2P\ 

L n r 2 J 



(8) 



where 2P is the effect of the imperfection of the ground plane which results from a finite value 
of the surface impedance Z. 

3. Evaluation of the Integral P 

The integral 

i/3\A)e- u(z + h) 



■P 



(u -f ifiA)u 
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J (\ P )d\ (9) 



must now be evaluated. Before proceeding, certain conditions and restrictions must be clearly 
stated. At a later stage, these may be relaxed somewhat. For the moment 

|A|<1, 0<arg A<J and ft>»l. 



Since 



J (x)=^H^(x)+H?\x)j 



and 

where flo 1} and H™ are Hankel functions of the first and second kind, respectively, it follows 
that 



> = *Af 8 



(u-\-ifHA)u 

Introducing the substitutions \=0 cos a and A=sin a, it is seen that 

i$A f * +i " cos ae-'W+v sln « H^jPp cos<*) 



m»(\ P )d\ (io) 



, = ^a r- 



sin ' cos ^ 



**. (11) 



Essentially P is now a spectrum of plane waves traveling away from the ground plane so that 
the vertical component of the propagation constant is ft sin a. Again, complex angles are 
included in the spectrum. When |/3p cos <*|>>1, the first term of the asymptotic expansion 
of the Hankel function can be employed. That is, 

H^iPp COS a)c±[——^ TV*P cos« 6 -*t/4. ( 12 ) 

\jirPp cos oj 
When this is inserted into eq (11), the equation for P becomes 

P=- e -*r,i 4 (JLY C' +t " (co*«)*e- iia "" 9 -"> da (13) 

2\2* P J J. tm sin ^ cQs ^ 

z z 

where R and are defined by z-\-h= R sin 6 and p= — i? cos 0. In the usual case, the separation 
p between transmitter and receiver is large compared to their respective heights, h and z, and 
consequently, 6 is slightly less than t. It can be seen that the exponential factor in the inte- 
grand of P is rapidly varying except for a region near <x=9^t. This is, of course, the saddle 
point of the integrand or the point of stationary phase. The important part of the integrand 
is in the region near the saddle point. In fact, this is the justification for employing the 
asymptotic expansion for the Hankel functions Ho l) ; the argument /3p cos a is always large 
in the region near the saddle point if 0p itself is large and z+h is less than p. 

The integral for P is now in a form where the saddle-point method of integration can 
be applied. The usual technique [23] is to deform the contour (which in the present case is 
along the negative imaginary axis, the real axis from to ir, and then along a line parallel to 
the positive imaginary axis) to a path of steepest descent. For example, if we let 

COS (a — d) = l— ir 2 
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and let r range from — <» to + » through real values, an integral of the type 

C G(cosa) e -^ Rcos(a - e) da (14) 

J-ioo 

where G (cos a) is slowly varying at a—0, is transformed to 

(2i)*6-** P " f (co8 a) e-v* T *dT. (15) 

In this deformation of the contour, account must be made for the singularities of the integrand 
that are crossed. In the case of poles, 2iri times the sum of the residues is added to the new 
integral. The final step in this classical procedure is to expand G(cos a)/Vl— ir 2 /2 in a power 
series in r enabling the integration to be carried out term by term. The leading term of this 
resulting asymptotic expansion is 

(2wi)KG(co* tiie-WtfR)-* (16) 

and succeeding terms contain (/3i?)~ % , (/3/?)~^, and so on. 

As is often the case, the contribution from the saddle point (which is the branch point in 
the original X plane) cannot be separated from the pole(s) of the integrand. For example, in 
the present problem, 

6r(cosa)=- -, (17) 

Sin COS — pr-^ 

has a pole at a=7r+a which is near the saddle point a— 6 since a is small and 6 is near 7r. In 
other words, the integrand is not slowly varying near the saddle point on account of the factor 
cos (a—a )/2. The integral is then of the form 



-s: 



TT + ico^-ipR cos(0-a) 

da (18) 



a—a 
COS ^ 



after the slowly varying factors have been separated out. The necessary modification of the 
saddle-point method to treat integrals of this type was devised by Van der Waerden [24] and 
Clemmow [25]. On making the usual deformation of the contour via the substitution cos 
(0—a) = l—iT 2 , the integral now becomes 

V " V 2 ^"^[cos^)] (19) 

It can readily be verified that no poles of the integrand are crossed in this deformation for the 
argument of a (or A) less than 45°. The integral to consider is now of the type 






2 +c 
where, for real r, x=fiR and 



-,r 2 

dr (20) 



Since f ^-^7=^) > {t fol l° ws that 
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Now it can be seen readily that 

i^ A =-(Sf Vxc (22) 

and an integration with respect to x from x to °° leads to 

e-*Vl=JJ(^yV~<fa, (23) 
or, after a change of variable, 

A=20/c)'V< : f _e~ z2 dz, (24) 

which holds for all c in the c plane with a cut along the negative real axis. This can be written 

A=^ e xc erfc^xc, (25) 

where erfc is the complement of the error integral given by 

erfcz =4= f e~ z2 dz y ( 26 ) 

where the integration in the z plane is directed from z towards the right to + °o . Therefore 

1= -i2Te- i P R e- w erfc(iw y *) ) (27) 



where 

and finally, since R=r 2 , 

where 
and 



w= 



-pRc=-2il3R |~cos (^2^)T> 

e- u ^ h) J (\p)d\ (28) 



(#XA) m _ uM} 



Jo (u +ipA) 
^i(irpy / £e~ w erfc(iw y2 )e~ i P r 2/r 2 , 



w ^( i+( Sr) 2 (29) 



i/3r 2 A2 , - ib 



While the derivation was carried out for the argument of A lying in the range to t/4 or b 
in the range to — tt/2, the formula is actually valid for all values of b if the above definition 
of the error function complement is used. The justification for this step is based on the principle 
of analytical continuation. If, in the initial derivation, the argument of A was allowed to exceed 
tt/4, we would find that a pole would be crossed in the deformation of the contour in the a 
plane. There would be a corresponding change in the integral A but, after a change of the sign 
of the variable r, the expression yielded for / would be unchanged in form if we are consistent 
in the definition of the complex error function. 

4. General Properties of the Solution 

The general behavior of the result will be considered now. For the moment the heights of 
the transmitter and receiver are assumed to be very small. That is (2+A)<<<p and, con- 
sequently, iv^p and r x ^r 2 ^p. The Hertz potential then has the form 

U = 2ce-^F(p), (30) 
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where 

F(p) = l-i (wp)Ke- p erfc (ip*)'. (3 1 ) 

The vertical electric field is then given by 

w l d dn (W 

E ^d P P d P (32) 

and, since /3p>>l and F(p) is a slowly varying function, it follows that 

S f = +2cP 2 e- i 0>p- l F(p). (33) 

For a highly conducting flat ground plane, [I#>|<<1], the electric field denoted by E? is 
given by 

ET = 2cfi 2 e- i ^p- 1 . (34) 

This result can be obtained by elementary methods, being just the free space radiation field of 
an electric dipole of moment 2c. The factor 2 arises from the imaging of the primary dipole in 
the perfectly reflecting surface. The ratio 

F( V ) = E 2 /E7 (35) 

is by definition the complex attenuation factor and is a measure of the efficiency of trans- 
mission of waves along the interface for a vertical-electric-dipole source. 

The behavior of F(p) is very dependent on the phase of p, and since 

P=M ,.__2p 4 .._*p (_*_)', (36) 

the surface impedance, Z, plays an important role. For example, in the case of a homogeneous 
ground of conductivity a and dielectric constant e, 

Z 1 A 1 V* 



A=- 



\120tt N 
with 



1 / IV* 



N jr *-i*/<* T. 



The refractive index, N, is usually large compared to unity with a phase near, but always less 
than —45°. Correspondingly, for a homogeneous half-space, the phase b lies between 0° and 
— 90°. This is precisely the case that has been treated in detail by Norton, who gives numerical 
results of F(p) as a function of \p\ for b varying from to — 180°. 4 In all cases for |p|>0, \F(p)\ 
was less than unity. In other words, the field E z was never greater than E™ for the perfectly 
conducting ground plane. For purposes of numerical work, the error-function complement 
can be expanded in an ascending power series, as is well known. This leads to the following 
series formula for F(p): 

F(p)=l-i(*py^e-*-2p+ { ^-^ ) +. . . , (38) 

which is absolutely convergent. 

The second term of this series, the quantity — i(irpY / -e~ p , is rather interesting if it is consid- 
ered by itself. The vertical electric field, E\ , which would result by ignoring the other terms, is 

Et= -2F/3 2 i(irpV)e- p (39) 

^A^-tU (40) 



« Norton used a time factor e- iut , whereas e iut is employed in the present paper. 
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This field has the characteristics of a Zenneck surface wave ; it varies inversely with the distance 
along with an exponential attenuation, and it has a phase velocity greater than that of light. 
By retracing the steps in the derivation it is not difficult to show that E\ is the residue at the 
pole a in the complex a plane. Taken by itself, E\ could easily exceed E™ by a great amount 
if b w r as near —90° corresponding to the phase of Z (or A) being near zero. This corresponds 
to a small conductivity in the lower half space. 

The above argument would indicate that the Zenneck surface wave is excited by a vertical 
dipole source on a homogeneous half-space. Of course, when the other terms in the series for 
F(p) are considered, the total field does not have the attenuation characteristics of a Zenneck 
surface wave. For example, if the asymptotic expansion for the complement of the error func- 
tion is employed, the following series for F(p) results: 

{P) 2p (2p) 2 {2pf (2pf "' K } 

which are also asymptotic and valid for large values of p and — 27r<6<0. This would indicate 
that the field E 2 at large ranges varies inversely as the square of distance, whereas the Zenneck 
surface wave varies inversely as the square root of the distance. This viewpoint is confirmed 
by the work of Wise [26] and Rice [27], who develop rigorous expansions for the Hertz potential 
on a homogeneous earth. 

The absence of any features associated with a trapped surface wave on a homogeneous 
earth is due to the restricted range of b (namely between and —t/2) for a homogeneous earth. 
It is of interest to note that Sommerfeld in his 1909 paper had an error in sign which was finally 
pointed out by Norton [3]. This error is equivalent to the interchanging of the b and —b. In 
other words, the complex conjugate of p is used in place of the proper quantity. Rolf [4] 
employed Sommerfeld's formulas along with the incorrect sign. Rolfs series expansions for 
F(p) would correspond to eq (38) with b replaced by —b. His calculations therefore correspond 
to values of F(p) for b from to +90°. The appearance of his curves is indeed weird. As 
will be indicated below, the function F(p) can exceed unity in certain instances when b is posi- 
tive. This is indicated in Rolf's curves. The irregularities and dips to zero, however, must 
be attributed to errors in computation. The series formula for F(p), although mathematically 
convergent, takes a great number of terms when p is greater than 5 or so. And, furthermore, 
the asymptotic expansion of F(p) is not accurate for the real part of p less than about 15. The 
anomalous features in Rolf's curves occur just where the accurate computation of F(p) would 
require tens or hundreds of terms in the power series formula. This viewpoint is confirmed 
when more elaborate calculations are carried out [21]. 

The nature of the attenuation function F(p) when b is positive can be best seen from an 
asymptotic development. This is obtained by noting that 

erfc(ip*)=>y= f e~ x2 dx=2—^= ( e~ x2 dx- (42) 

The latter integral asymptotically vanishes for p^>co if the real part of ip 1/2 is less than zero. 
Therefore 

for 27r>&>0 and|^|>>l. Utilizing this result, it then follows that 

F(p) = l—i^/irpe~ p erfc(ip^) 

^—2i^irpe~ p 



(44) 
1 1-3 1-3-5 



2p (2p) 2 (2pf ' ' • 
in the asymptotic sense. The term —2i^wp e~ p has all the characteristics of a surface wave. 
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It is not present in the asymptotic development for F(p) when b is negative. Of course, for 
dissipative media such as a stratified conductor, the real part of p is large and the surface-wave 
term almost vanishes compared to the term — l/2p. 

The most interesting case is when the surface impedance is purely inductive. For example, 
if a flat perfectly conducting ground plane is coated with a thin layer (of thickness h) of pure 
dielectric (with refractive index JV), then the surface impedance, Z, is given to a good approxima- 
tion by [19] 

Z^Wfc[l~]. (45) 

Another example is when a flat metal surface is corrugated in the form of periodic rectangular 
slots of width d and depth D with a separation w between slot centers. The surface impedance 
for the magnetic vector parallel to the directions of the corrugations is then 

Z=i-ifotan/3Z?, (46) 

provided that there are at least five corrugations per wavelength along the surface [13]. For 
D<X/4, the surface impedance is purely inductive, and since 

p=-^(Z/, ) 2 , (47) 

p+i|p| and b=ir/2. The attenuation function then lias the asymptotic development, 

F( P ^-2e^me-^+^+^-^- • • , (48) 



where 



b|=UandX=|Z|. 



The vertical electric field at z=0 is then given by 

E z ^2cP 2 e-^p- l F{ V ). (49) 

The surface wave term now varies as 

which is not a Zenneck surface wave since the phase velocity, v, given by 



_ 1 1 

V 1+^/(2^) V^' 



(50) 



is less than the velocity of light, l/vWo. 

If the surface is purely capacitive such that Z=—iX', the argument of the error function is 

i p V=e-«'*y[M2(X f /r, r, (51) 

so that erfc(ip^) can be expanded directly to lead to the asymptotic representation 

F ^-2M+2||p + # ' ' ' (52) 

which contains no surface-wave term. 
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5. Height-Gain Functions 

It is desirable now to return to the general formula for F(p) which is a (modified) saddle 
point approximation for the field of a vertical dipole over the impedive surface. For h = 0, it is 

W=l-i^e- w erfc(iw^), (53) 

where 

w= ^( 1+ i 2 ) 2 ^( 1 +4) 2 ' 

For small heights and receiver it is of interest to develop Win a power series in z. For example, 

TO=mO)+<f )+| (*£);_+• • • • (54) 



Now 



iind 



dW T p../ — „ . ,. x/ ~\dw . rN 



Therefore 



and, to a first order in 



dz P \^ApjAp 



f^ =ikAW(0) = ikAF(p) 



W(z)o*F(p) [1+iifcAs]. (56) 



The generation of higher-order terms using this method is cumbersome, and it is really better 
to return to the original integral for the Hertz potential n(s) given by 

U(z)=2cr X ^^e-^d\ (57) 

Jo u + ipA 

for h=0. Expanding the exponential, this is rewritten 

Jl(z)=2j"ll-uz+uV/2- . . . ]£$j£dk (58) 

=2c[n ( 0) + ^n ( 0) + ^(^)^+...} 

Now since II (z) satisfies the wave equation, tf/bz 2 can be replaced by 

1 b b 



32 — *- p sr 

p Op Op 



so that 



n (z )=[i + ^A 2 -^+I|p|)+...]n ( o). 



(59) 



Since 

ll(0)=2ce- i ^p- l F(p), (60) 
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the contribution from the second-order term is 

where the prime(s) indicate a derivative(s) with respect to p. It is not difficult to see that 
for 2w <^6<C0, F(p) is slowly varying compared to other factors and the term in curly brackets 
is much less than unity if /3p>>l and p 2 >>0 2 . In the case of the inductive surface (i. e., 
27r>6>0), it is convenient to write 

F(p)=F s (p) + F a (p), (62) 

where asymptotically 

F,(p)&-2iifipe-* (63) 

is the surface wave term and 

*fr>"-^i|dp- • • • (64) 

is the remaining (nonsurface) contribution. [For the case — 2.t<&<0, F(p)^F a (p)]. 
Now 



Fa(T) 

and 



>@ 



but 

F.{p)~ 2 
and 

TO ( jfi**\ 2 

F,{p)~\ 2 A 

Therefore, the Hertz potential H.„(z), corresponding to the surface-wave field, is given by 

n,( 2 )sn,(o) [i+i/3A2+ <^M!+ • • •] (65) 

s n s (0)e^ 
whereas 

n o (2)sn a (0)[l-f-ii8A2]. (66) 

In the case of a purely inductive surface (Z=iX), the exponential factor for the surface 
term is given by 

e i/3A« =g -(Z/5/, )z = g-« M X« ) (67) 

indicating a rapid decay of the field with height. 

By a slight extension of the above, the case of both transmitter and receiver raised yields 

U s (z,h)^U s (0,0)e i ^ z+ " ) (68) 

and 

n a (s,A)sn.(0,0)(l+ip , A3)(l+i/3AA). (69) 
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Still another approach to the effect of raised antennas, which is particularly appropriate 
for large heights, is now considered. In general, 

==- +- 2P, (70) 

c n r 2 

where 

P=l [l-F(w)] — (71) 

w L r 2 

with F(w)=F a (w) + eF s (w), where e=0 for -2tt<6<0, e=l for 2tt>6>0, and w=p(l + C/A) 2 , 
0= (z-\-k)/r 2 . Asymptotically, 

F s (w)=—2i-y/rw e~ w (72) 

and 

F a (w)=-~-^- 2 --^ (73) 

aK J 2w 4w 2 8w* v ' 

After some algebraic manipulations, 

n=n a + € n s 

Ua^e-Wi . C-Ar^ J 1 1.3 1-3-5 *\ e'W* (7A) 

c^'fi ^O+A r 2 + \ p{l + C/A)* 2p\l + C/Ay ±p 3 (l + C/A) 7 ' * "/ r 2 { ] 

and 

S-^^V^^^]^- 2 - (75) 

c A+C r 2 

The first term in the expression for IT a is, of course, the direct or primary influence. The 
second term is the secondary influence modified by a Fresnel reflection coefficient. The terms 
containing inverse powers of p can be considered as corrections to geometrical optics. The 
surface-wave term is absent when b is negative. When b is positive, it is a finite contribution 
to the total field IT, and for a highly inductive surface (i. e., 6 = 90°) it is predominant. 

The vertical electric field is simply given by 

E z ^p 2 S 2 Il=p 2 (l-C 2 )n for fir x and £r 2 >>l, 

and furthermore, r 2 ^ri—2hC. 

6. Conclusion 

It has been shown that the nature of the field radiated from a dipole over a flat surface 
depends markedly on the complex value of the surface impedance Z. For example, when the 
phase of Z is between and 45°, the derived expressions for the field correspond to those of 
Norton and others for a dipole over a homogeneous flat ground. When the phase of Z exceeds 
45° the radiated field becomes trapped by the surface. This effect is most pronounced when 
Z is purely imaginary, corresponding to an inductive surface. A similar unifying treatment 
for a dipole over a curved spherical surface will be given in a sequel to this paper. 



The author thanks Kenneth A. Norton for his very helpful suggestions and Howard 
Bussey for his constructive criticism. 
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